An analytical approach to sorting in periodic potentials 
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There has been a recent revolution in the ability to manipulate micrometer-sized objects on 
surfaces patterned by traps or obstacles of controllable configurations and shapes. One application 
of this technology is to separate particles driven across such a surface by an external force according 
to some particle characteristic such as size or index of refraction. The surface features cause the 
trajectories of particles driven across the surface to deviate from the direction of the force by an 
amount that depends on the particular characteristic, thus leading to sorting. While models of this 
behavior have provided a good understanding of these observations, the solutions have so far been 
primarily numerical. In this paper we provide analytic predictions for the dependence of the angle 
between the direction of motion and the external force on a number of model parameters for periodic 
as well as random surfaces. We test these predictions against exact numerical simulations. 

PACS numbers: 05.60.Cd,66.30.-h,82.70.Dd,05.40.-a 
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I. INTRODUCTION 

Transport of particles driven by external forces across 
modulated potential surfaces, and the associated ability 
to sort mixtures of particles according to some feature 
sensitive to this modulation, has attracted considerable 
recent interest 

[EISIIillSllGimill The 

surface modulation may be periodic or random, and it 
may consist of traps, obstacles, or a combination of both. 
Typical sorting parameters include particle size, particle 
index of refraction, and particle mass. The underlying 
mechanism in all cases relies on the fact that while the 
particles attempt to follow the external force, the traps 
or obstacles induce systematic deviations in the trajec- 
tories so that there is a non-zero average angle between 
the trajectory and the external force. If this non-zero 
angle depends on some characteristic of the particle such 
as size, then particles that differ in this characteristic 
emerge at different angles and can therefore be sorted in 
a very effective way. 

There are a number of experimental demonstrations 
of these capabilities P, 0, 0, IE 1 and there are also 
fairly extensive numerical simulations of various model 
formulations 0, 0, 0, ITTl IT^ | that are quite successful 
in explaining the experimental phenomena. These sim- 
ulations to some extent clarify the roles of a number of 
physical variables such as, for example, the direction of 
the external force, surface geometry, and ambient tem- 
perature, that affect the experimental outcomes. On the 
other hand, analytic results are relatively rare and 



yet it is undoubtedly useful to have predictive analytic 
formulas, especially with a view to optimizing the sorting 
mechanism. In this paper we take a step in this direction, 
deriving approximate results that are shown to reproduce 
some of the important earlier numerical results. While 
our results are limited to certain parameter regimes, we 
believe that they may prove useful in designing experi- 
mental sorting potentials and in identifying likely regions 
of parameter space for more detailed numerical study. 



In Sec. |n] we describe the model and define the quan- 
tities to be calculated, specifically, the angle between the 
average velocity of the particles and the external force. 
In our discussions we include periodic potentials as well 
as random potentials. In Sec. IIIII (accompanied by an 
Appendix) we present approximate equations for the av- 
erage velocity valid for sufficiently strong forces and/or 
high temperatures, and discuss how these approxima- 
tions might be continued to further orders than those re- 
tained here. In Sec. lIVI we compare our theoretical results 
with those of numerical simulations for periodic poten- 
tials, and in the course of this comparison we refine our 
approximations so as to eliminate an unphysical outcome. 
With this adjustment we improve the theoretical predic- 
tions. In Sec. we carry out a similar comparison for 
random potentials, and again find very good agreement 
between theory and numerical simulations. We summa- 
rize our results and indicate some future directions in 
Sec. EH 
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II. THE MODEL 

We consider identical non-interacting particles moving 
across a surface described by a two-dimensional potential 
V(x, y) of unit height or depth and unit period in both 
directions. The equation of motion for the particles is 
given by a Langevin equation for each component of the 
particle displacement 
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d_ 

dx 
d_ 
dy 



V(x,y)+Fcos6 + Vzf( x (t) 
V(x,y) + Fsm6 + Vzf£y(t). 



(1) 



The dots denote time derivatives, T is the dimensionless 
temperature, and the thermal fluctuation terms are 
Gaussian and 5-correlated, 



(O) =SijS(t-t'). 



(2) 



Note that the equations are in dimensionless form, with 
the scaled temperature T and force magnitude F given 
in terms of physical parameters as in Eq. (4) of [Icj . 
The system is assumed to be overdamped so that inertial 
terms x and y have been dropped from the equations of 
motion. The constant external force vector is 



F = Fcos#i + Fsin#j, 



(3) 



such an approximate solution in the next section, and 
also generalize the results to the case of a random poten- 
tial by averaging over the disorder to obtain (v) . 

Of particular interest in sorting applications is the de- 
flection angle a of the velocity from the external force 
direction, given by the relation [T(| 

(«!_) ~ (vx) sin 9 + (v y ) cos 9 
tana = -, — r = — ; — ; ; — — — — . (7) 

\v x ) cos 9 + {Vy) sm9 

(If the potential is random, each velocity component 
{v±), (v x ) etc. in this formula should be replaced by 
its average over the disorder: (v±), (v x ) etc.). We will 
explore the behavior of this angle through that of each 
component (v x ) and (v v ) of the average velocity, and ex- 
amine the accuracy of our approximations in capturing 
these behaviors. 

To avoid overly complicated formulas, we adopt the 
convention throughout the remainder of this paper (ex- 
cept in the Appendix) that the notation v (or its compo- 
nents v x , v y , v±, vn) stands for the averaged quantity (v) 
(or, respectively, (v x ), (v y ), (v±), (v\\))- If the potential 
is random, v (or its components) stands for (v) (or its 
respective averaged components). Averages of all other 
quantities will be denoted explicitly using angle brackets 
and overbars as appropriate. 



and all particles are assumed to start at the origin. 

Our goal is to calculate the direction and magnitude 
of the velocity of the particles averaged over the ther- 
mal fluctuations. The long-time limit of this velocity is 
denoted by (v), and is decomposed into Cartesian com- 
ponents as 

(■v) = (v s )i+(v v )j. (4) 

The angle brackets denote averaging over the thermal 
noise. If the potential V(x, y) is random, then we also 
perform an average over realizations of the potential, de- 
noting such an ensemble average by an overbar: (v). To 
calculate the average velocity for a periodic potential, we 
write the Fokker-Planck equation for the concentration 
c(x, t) of the particles obeying the equations of motion 

09. 



| + V.[(F 



VV) c]-TV 2 c = 0, 



(5) 



with initial condition c(x, 0) = <5(x). If this equation 
could be solved for the concentration c(x, t) at time t, 
then the desired average velocity vector follows from 



lim 



fix 
~dt 



= lim — 

t— >oo at 



J dxxc(x, t). 



(6) 



While an exact solution of Eq. © is not in general pos- 
sible, an approximate solution for the concentration may 
be derived using a number of methods. We implement 



III. APPROXIMATE SOLUTIONS 

The experimentally interesting regimes involve an ex- 
ternal force of magnitude sufficiently larger than the well 
depth (F > 1) so that particles do not easily become 
trapped. If a particle did become trapped, it would have 
to be extracted by a sufficiently large thermal fluctuation 
to continue moving across the surface. The particles of 
interest in a sorting context are those that do not be- 
come trapped during the course of the experiment. On 
the other hand, if the force is too large, it simply drags 
the particles along, the features of the potential become 
essentially invisible, and particle separation, which relies 
on the effect of the potential on the particle trajectories, 
does not occur. At the same time, the temperature of the 
system should not be so high as to obliterate the features 
of the potential. A temperature that is too high would 
again lead to an uninteresting situation in that particle 
separation would again not be observed. The regime of 
interest is thus that of an external force that is large but 
not too large, and a temperature that is as low as possible 
in theory and in practice. 



A. First order approximation for periodic 
potentials 

In the Appendix we show how a systematic large 
force and/or high temperature approximation may be 
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obtained. The first-order approximation leads to the av- 
erage velocity vector 



v = F 



k 2 k 



Tk 2 - ik ■ F 



Q(k). 



(8) 



where for periodic potentials the function Q(k) is defined 
by 



Q(k) = V(k)V(-k). 



(9) 



Here fc 2 = k.k and V is the Fourier transform of the 
potential V(x) as defined in Eq. (|A1() . 

For later comparisons, it is useful to exhibit the result 
of this approximation for the simplest one-dimensional 
version of our problem, one with a simple cosinusoidal 
potential. This case has been studied extensively (see, for 
example, Chapter 11 of ,13]), and exact solutions may be 
found for the average velocity. The equation of motion 
is 

dx i 

— = C7 + 27rsin(27ra:) + V2T£(t). (10) 

We call the forcing U instead of F because for the sim- 
plest separable potential in two dimensions we will be 
able to identify the average speed v(U) = (x) with the 
components of the two-dimensional velocity via the re- 
lations v x = v(Fcos8) and v y = v(FsinO), thus pro- 
viding a way to determine the deflection angle in a two- 
dimensional case via a one-dimensional calculation. The 
result is 



v(U) « U 



2ir 2 U 



4tt 2 T 2 + U 2 ' 



(11) 



A number of points are noteworthy. Firstly, a simple 
sinusoidal potential has minima (which we think of as 
traps) as well as maxima (which we think of as barriers) 
relative to a flat landscape 0. Secondly, we observe 
that the net effect of the potential in one dimension is to 
slow down the particles. Thirdly, when the associations 
Vx = v(FcosO) and v y — v(Fsm9) are appropriate, we 
can immediately see that the deflection a is associated 
with the fact that the "slowing down" is different for each 
component for most angles 9. 

The technique yielding the approximate formula JSJ in 
the case of a periodic potential may readily be extended 
to random potentials, as simulated numerically in |14| . 
for instance. Indeed, the method of deriving Eq. (JSJ) was 
originally developed for the calculation of particle con- 
centrations when advected by potentials which are ran- 
dom in space and also in time 0, 0] . The possibil- 
ity of sorting on random potentials can thus be examined 
using analytic approximations. 



where E(k) is the "energy spectrum" of the potential. 
The energy spectrum is defined as the Fourier transform 
of the (disorder-averaged) correlation function of the po- 
tential V (x') V (x' + x). The disorder is assumed to be 
homogeneous, so that the correlation function depends 
only on the difference vector x; recall that the overbar 
denotes averaging over the ensemble of random poten- 
tials. 

It is again useful to exhibit the first order approxi- 
mation for the simplest one-dimensional version of the 
problem. The equation of motion is 



dx 
~dt 



u 



dx 



V(x) 



(13) 



where V(x) is now a zero-mean random modulation with 
spectrum E(k). The one-dimensional velocity averaged 
over realizations of the noise and of the random po- 
tential is then approximately given by 



v(U) = U- 



u_ 

2^ 



dk 



k 2 E{k) 
T 2 k 2 + U 2 ' 



(14) 



A more explicit result requires the specification of the 
correlation function. 

It is an interesting general result that regardless of 
the specific form of the correlation function, the aver- 
age deflection angle is zero when the random potential in 
our two-dimensional scenario is isotropic, that is, when 
the correlation function V (x') V (x' + x) depends only 
on the magnitude |x| of the difference vector and not on 
its direction (see, for example, Eq. (14) of 01). In this 
case the energy spectrum E(k) depends only on the mag- 
nitude k — \k\ of its argument. Using l|12(l . we consider 
the wave-vector integral in (JHJ by writing k in terms of 
a component k\\ = kcoscf) parallel to the force F, and a 
component fej_ = ksm(f) perpendicular to F, where <p is 
the angle between k and F. With the two-dimensional 
integral written in polar coordinates, the average perpen- 
dicular velocity is found to be 



Vj_ 



(27T) 2 



dk k 



d</> 



k 3 sin ( 



Tk 2 — ikF cos i 



■E(k). 



The isotropy of the potential means that the energy spec- 
trum is independent of the angle <f>, and therefore the 
integration over cf> may be performed in l|15f) . giving the 
result v± = 0. We conclude that for isotropic random 
potentials the deflection angle is zero. This means that 
sorting by deflection angle is not possible in isotropic ran- 
dom potentials. 



B. First order approximation for random potentials 

The appropriate kernel in formula © for a random 
potential is 

Q(k) = 4^ 2 £(k) (12) 



C. Higher order approximations 

It is useful to explore the behavior and magnitude of 
the second and higher order approximations for the ve- 
locity. We do not do this in full generality, but only for 
the simplest one-dimensional problems (|10fl and (|13|) . 
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Following the procedure outlined in the Appendix, a 
second iteration is somewhat cumbersome but straight- 
forward, and leads in the periodic case to the approxi- 
mate velocity 



v(U) w U 



2tt 2 U 



-2tt 4 



47r 2 r 2 + u 2 

U 3 - 20n 2 UT 2 
(U 2 + 4ir 2 T 2 ) 2 (U 2 + 16tt 2 T 2 ) 



(16) 



The last term is the second order correction. In a sub- 
sequent section, where we make comparisons with exact 
results, we comment on some higher order terms in this 
series. Here we simply point out that for T = the 
full infinite series sums to the known exact solution Hal 



v(U)/U= y/1 - 4tt 2 /[/ 2 . 

In the case of a one-dimensional random potential, con- 
tinuing the iteration to the next order extends Eq. (f 1 4|) 
to 



v(U) = Ti- 



ll 
2^ 



dk 



k 2 E(k) 



1 



(27T) 2 



— oo 

OO /'OO 



T 2 k 2 + U 2 



oo J —oo 



dp dqp 2 E(p)q 2 E(q) 



(iU - Tp)(iU - T{p + q))(iU - Tp) 
i 



{iU - Tp)(iU - T{p + q)){iU - Tq) 



(17) 



Finally, while we recognize that the approximation 
methods based on the truncation of a series precludes a 
straightforward determination of the regimes of validity, 
one might expect that high temperature and/or strong 
external forcing (T ^> 1 and/or F ^S> 1) would be suffi- 
cient to make the approximations presented in this sec- 
tion reasonably accurate. As we shall see later, such an 
assertion requires some caveats, but numerical simulation 
results help resolve the issues which arise. 



IV. RESULTS FOR PERIODIC POTENTIALS 

In this section we test the approximate results against 
exact simulation results for two-dimensional periodic po- 
tentials, explain why large F is not necessarily sufficient 
for agreement at the orders developed in the last section, 
and modify our approximations nonperturbatively so as 
to correct for the source of disagreement. 

If the potential is 1-periodic and even in both x and 
y, such as the one used in [ljj, it may be expanded in a 
Fourier series as 



M M 
n— m— 



, cos(27rnx) cos(27rm?/), (18) 



where M may be infinite, although in our calculations we 
retain only a finite number of modes in each direction. 



After some algebra, Eq. JSJ) yields explicit expressions 
for the average velocity in the x and y directions as sums 
over Fourier modes, 

2 M M 

v x = Fcos9 — 2J a nmd n d m Fn(n 2 + m 2 ) 

n— m— 

n cos 9 — m sin 6 
Att 2 T 2 (n 2 + m 2 ) 2 + F 2 (n cos 9 — m sin 9) 2 
n cos 9 + m sin 

+ 



47r 2 T 2 (n 2 + m 2 ) 2 + F 2 (n cos 9 + msin#) 2 

2 M M 

F sin 9 - — ^2 X! a nrn d nd m Fm(n 2 + m 2 ) 

n— m— 

77i sin 9 — n cos 9 



+ 



47r 2 T 2 (77 2 + to 2 ) 2 + F 2 (ti cos 9 — msiiit 

m sin 9 + 77 COS 9 



AiT 2 T 2 (n 2 + m 2 ) 2 + F 2 (ncos9 + msin( 



(19) 



where the coefficients a nm depend on the specific poten- 
tial, and 



d n = 1 + 5 n o 



2 if 77 = 
1 if 77 > 0. 



(20) 



Note that to this order of approximation the Fourier co- 
efficients a nm appear only as a 2 , so that the resulting ve- 
locities are the same regardless of the overall sign of the 
potential, e.g., whether the potential consists of traps or 
of obstacles. 

To examine the usefulness of the approximate formulas 
(|19fl . we choose the periodic potential used in [Io| . 



V(x,y) 



1 



1 



(21) 



with g(x,y) = 5 [cos(27rx) + cos(27ry) — 2B], and with 
various values of the parameter B (associated with dif- 
ferent particle sizes jldj). The negative sign in the nu- 
merator indicates that we are considering traps (but, as 
noted above, this sign does not matter for the assessment 
of the validity of the approximation) . For each value of 
B, the potential is first expressed in terms of its Fourier 
series (IT%|) to determine the Fourier coefficients a nm . 

Figure n shows the angle tana defined in Eq. (JJJ plot- 
ted against tan# for F = 8 and T = 0.1, and with 
B = 0.9. Note the degree of agreement between the exact 
and approximate results, especially the positive and neg- 
ative regions of deflection angle which reflect the "terrace 
phenomenon" or "locked-in" states observed in experi- 
ments 0, U S S S 13 • Note also the accord in the order 
of magnitude of the deflection. These indications point to 
the qualitatively satisfactory performance of the approx- 
imate formula 10 even at the low temperature T = 0.1. 
Also of interest here is the dependence of the results on 
the number M of Fourier modes retained in the expan- 
sion l|18|) of the potential. The shape of the potential 
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-0.020 1 ' 1 ' 1 ' 1 ' 1 ' 1 

0.0 0.2 0.4 0.6 0.8 1.0 

tan 6 

FIG. 1: The dependence of the deflection angle a on the 
forcing angle 6 at temperature T = 0.1, and for parameters 
F — 8 and B — 0.9. Numerical simulation results (solid 
diamonds) agree well with theoretical predictions using M — 
6 modes (solid line) in the Fourier series 1181 . For comparison, 
the theoretical prediction using only M = 1 modes is also 
shown (dashed line). 

(|21[l is well approximated with M = 6; however our re- 
sults with M as low as 1 show that the correct order of 
magnitude of the deflection can be predicted using only 
the lowest Fourier harmonics of the potential. 

The Fourier decomposition provides an opportunity to 
understand the role of various symmetry contributions 
to the sorting capability of the surface. For this purpose, 
we examine two special cases with M = 1. The potential 

V a (x, y) = cos(2ttx) + cos(27ry) (22) 

corresponds to having a 01 = am = 1, with all other a nm 
being zero. This separable potential was used in ^1 an d 
recently in 0], and we discuss it further subsequently. 
For the parameters used in Fig. ^ the potential V a yields 
deflection angles which are negative for all forcing direc- 
tions. In contrast, the potential 

Vb{x, y) — cos(27rx) cos(27ry), (23) 

corresponding to an = 1, with a nm = otherwise, gives 
positive deflection angles for all forcing directions. It is 
the (appropriately weighted) combination of the poten- 
tials V a and Vb in the Fourier series of 1)21(1 which gen- 
erates the crossover from negative to positive deflection 
angle in the M = 1 curve of Fig. ^ 

In Fig. we show results of numerical simulations of 
sorting in the separable potential V a at temperature T = 
0.1 and for various magnitudes of the forcing F. As noted 
following Eq. 1)10(1. separable potentials allow the appli- 
cation of one-dimensional results to the two-dimensional 
case by setting v x = v(Fcos8), v y = v(Fsm9). The 



dashed curves in Fig. [S] show the deflection angle pre- 
dicted using the first order approximation lfTT|) . While 
the agreement is quite good at the higher angles 8, it is 
clearly not quantitatively satisfactory at low angles. The 
second order approximation from Eq. 1(16(1 yields the solid 
curves, and further improves the agreement at the higher 
values of 9 and pushes this agreement toward lower 9, but 
does not greatly improve the low-# situation. The diffi- 
culty can easily be traced in the case of the separable 
potential and, by inference, for more complex potentials, 
by recalling that v y = v(Fs'm9). Clearly, at low angles 
the argument of v y is necessarily small no matter the 
magnitude of F. Equations (|TTJl and (|T7jl make it clear 
that at sufficiently small argument the one-dimensional 
approximate velocity becomes negative. This unphysical 
result is reflected in the misbehavior of v y at small 9 and 
leads to the disagreements observed at small angles. 

To test this hypothesis, we introduce an enhanced 
approximation, which consists of calculating v x and v y 
from the first-order or second-order iteration approxi- 
mation, as we have done so far, but then replacing v y 
by max[iij,,0] and v x by max^^O]. This eliminates the 
unphysical negative values of the components, replacing 
them by zero, without affecting positive outcomes. We 
shall call this the adjusted truncation, and its effect on 
the predictions is to cut off the curves to the left of the 
dotted line tana = — tan 9 in Fig. [5] The agreement 
between the adjusted truncation and the numerical re- 
sults is clearly very good , with the remaining differences 
arising from the abrupt replacement of negative velocities 
by zero velocities in the adjusted truncation (and, for the 
second order case, a very small region near 9 — 0, where 
a positive deflection angle is erroneously predicted) . 

We end this section by repeating our earlier assertion 
that the perturbation theory becomes more accurate with 
increasing temperature. In fact, the series expansion de- 
veloped in the Appendix is a perturbation series about 
the limit of infinite forcing, and non-zero temperature 
acts to regularize this series at finite forcing. The ac- 
curacy of truncations of the series thus improves with 
increasing force and with increasing temperature. Al- 
though the sorting capability of the system decreases 
with increasing temperature (and with increasing force), 
it is useful to display some higher temperature results 
explicitly. We do so in the one-dimensional case. By 
continuing the iteration method introduced in the Ap- 
pendix, we write successive approximations to the aver- 
age velocity as 



N 

v ( u ) = U + 2nJ2 v n, (24) 

n=l 
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FIG. 2: Deflection angle for separable two-dimensional po- 
tential 122H with T = 0.1 and F = 15, 25, 50 from greater to 
smaller deflection angles. Dashed curves: first order approx- 
imation 111H ; solid curves: second order approximation 1161 1: 
dotted line: cut off for approximation curves in the adjusted 
truncation; symbols: numerical simulation results. 



with 



u 



vi = 



V 2 



2(C/ 2 



+ r 2 ) 

U 3 - 5UT 2 



v 3 



8(U 2 + T 2 ) 2 (U 2 +4T 2 ) 

U(U 4 - 24U 2 T 2 + 23T 4 ) 
16(U 2 + T 2 ) 3 (t/ 4 + 13U 2 T 2 + 36T 4 ) 



(25) 



Here we write U = U/(2tt) to keep the formulas simple; 
note that v\ and v 2 have already appeared in Eq. (|16|l . 
Figure |21 show v/U vs T at forcing value U = 3 as the 
number of terms retained in the series is increased. Note 
that at high temperatures, even the simple first order 
(N = 1) truncation already gives accurate results. At 
lower temperatures, retaining higher order terms in the 
approximations gives results which are closer to the ex- 
act values. As the temperature approaches zero for this 
U value, the approximations all fail; this effect depends 
on the forcing U, and at higher values of U all the ap- 
proximations remain accurate even at T = 0, as noted 
following Eq. ((Tfj|l . 



FIG. 3: Results for v/U for the one- dimensional potential 
with forcing U — 3, plotted as a function of temperature T. 
Symbols show exact values; the approximation curves use v 1 -^ 
for v, as denned in equation 12 U . with coding: TV = 1 (solid), 
N = 2 (dashed), N = 3 (dotted), N = 4 (dash-dotted). 



Recall that sorting is not possible in an isotropic ran- 
dom potential. We must therefore choose an anisotropic 
potential. In particular, we consider a two-dimensional 
(separable) potential generated by adding together inde- 
pendent, zero- mean modulations in the x and y direc- 
tions, 



V(x,y) = V 1 (x) + V 2 (y), 
with correlation functions given by 



(26) 



Vi{x')Vi{x' + x) = Ei{x), 
V 2 (y')V 2 (y' + y) = E 2 (y), 



W)W) = 0. (27) 
The energy spectrum of this potential has the form 

E(k) = 2n6(k y )E 1 (k x ) + 2TTS(k x )E 2 (k y ), (28) 

where E\ and E 2 are the one-dimensional spectra given 
by the Fourier transforms of V\ and V 2 . Using spectrum 
(|28|l in the formula © yields the average velocity com- 
ponents 



V. RESULTS FOR RANDOM POTENTIALS 

We next test our approach for random potentials, for 
which the first order approximation to the velocity is 
given in general in Eq. © and Eq. 1)12(1 . and in one di- 
mension in Eq. I|14|l . The second order iteration for a 
one-dimensional random potential is shown in Eq. 81711 . 



Fcos6 



Fsin( 



1 

2^ 



elk 



k^Exi^FcosO 
T 2 k 4 + k 2 F 2 cos 2 6>' 

k 4 E 2 (k)F sin 6 
2W__ T 2 k 4 + k 2 F 2 sin 2 e»' 



— OO 

oo 



(29) 



We note that these formulas imply that the average veloc- 
ity component perpendicular to F is non-zero in general 
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-0.25 1 1 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

tan 6 



FIG. 4: Deflection angles in the two-dimensional separable 
random potential with Gaussian energy spectrum in each di- 
rection (s = 5, A = 4), temperature T = 0.2 and external 
forcing F = 1 (lower results) and F — 2 (upper results). The 
dashed curves show the first-order approximation, and the 
solid curves the second-order approximation. The adjusted 
truncation cuts off both approximations at the dotted line. 
The numerical simulation results are shown as symbols. 



(even if E\ = E 2 )' 



VI. CONCLUSION 

In summary, we have derived approximate analytic or 
quadrature forms for the average velocity of particles 
moving in a periodic potential or random potential. The 
average velocity in general deviates from that of the ap- 
plied constant external force, the deviation depending on 
some particle characteristic such as size. This can then 
be used to sort particles that differ in this characteristic. 
A systematic perturbation series valid for large external 
forces and/or high temperatures is shown to capture the 
behavior observed in experiments and numerical simula- 
tions in physically interesting parameter regimes when 
the angle between the external force and the crystallo- 
graphic x axis (or y axis) is relatively large, e.g., when 
the force lies near the crystal diagonal, even when the 
perturbation series is truncated at low orders. The re- 
sults are not nearly as good when the force lies near one 
of the crystal axes, because the truncated perturbation 
series can then lead to unphysical negative velocity com- 
ponents. We have proposed an adjustment to the simple 
perturbation expansion whereby negative velocity com- 
ponents are set to zero, and have shown that this adjusted 
truncation scheme leads to very good agreement with nu- 
merical simulation results even for low temperatures. 

Further directions in this work are plentiful, and here 
we list just a few. A first direction would be a gener- 
alization to reduced-symmetry potentials. For example, 
in 01 we considered a potential of the form 



-v x smt 



F cos 9 sin ( 



2tt 
E 2 (k) 



dk k 4 



Ei(k) 



T 2 k 4 + k 2 F 2 sin 2 



T 2 k 4 + k 2 F 2 cos 2 6 
(30) 



The separable potential (|2(j|) is special because it reduces 
the two-dimensional sorting problem to motion in un- 
coupled one-dimensional potentials in the x and y di- 
rections. The results <|29|) for the two-dimensional case 
follow from Eq. (|14|l by noting that v x = v(Fcos0) and 
v y = v(Fsin6). 

A comparison between the first-order adjusted trun- 
cation results and numerical simulations for a two- 
dimensional separable potential with correlation func- 
tions 



/2ny 

V(x, y) — cos(27rx) cos I — — 
V V 2 



(32) 



which has different length scales in the x and y direc- 
tions. Another extension would be to the approximate 
calculation of the diffusion tensor 



lim — ((a;, - (x,)) (x 3 



(xj))) , (33) 



(see, e.g., Eq. (4.24) of |l5j). in particular to compare val- 
ues for diffusion in the transverse and parallel directions 
to the direction of transport [l^] ■ Reimann et al. sup- 
ply an exact expression for the effective diffusion in one- 
dimensional problems, which should also be applicable 
to the two-dimensional case with a separable potential. 
These and other related investigations are in progress. 



/ 27rA 



D -x 2 /2A 2 



My) = 



V/2A 2 



(31) 

is presented in Fig. While the agreement is not as 
good as in the periodic case, the theory clearly captures 
the numerical behavior rather well, and the effect of the 
order of the approximation is again manifest. 
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APPENDIX A: ITERATION 

In this Appendix we outline a systematic approxima- 
tion scheme that leads to our main results in Sec. IIIII 
starting from the Fokker- Planck equation JSJ. We call 
u(x) = — VV(x) for convenience. The case where u is 
a Gaussian random field perturbing the strong external 
bias F has been examined in, for example, Sec. 4.2.2 
of 0| . In our case the field u may be random or periodic, 
and we mimic the derivation of the systematic expansion 
of ^3 to determine our results for the long-time average 
velocity. 

Defining the (spatial) Fourier transform c of the con- 
centration field c by 



c(M) 



rfxe- ,kx c(x,«), 



(Al) 



the Fokker-Planck equation is transformed to the equa- 
tion 

dc % f 

— +ik-Fc+-^-^ J dpk-u(p)c(k-p,t)+k 2 Tc(k,t) = 0, 

(A2) 

with initial condition c(k, 0) = 1. Note that the Fourier 
transform of u is 



u(p) = -ipV(p). 
As a consequence we have 

u(0) - 0, 



(A3) 



(A4) 



which we will use below. 

We also perform a Laplace transform in time, defining 



/>oo 

c(k,s) = / dte- st c(k,t). 
Jo 



(A5) 



After Laplace transforming Eq. i|A2(l and using the initial 
condition, our goal becomes the solution of the integral 
equation 



c(k, s) =P a (k) 



(2tt)2 



P s (k) / dpk-u(p)c(k-p, S ), 

(A6) 

where the propagator (or "free Green function" |15|) is 
defined as 



^(k) = 



1 



s + k 2 T + ik-F' 



(A7) 



The following limiting behavior of the propagator will be 
important later: 



limsPfk)-! ifk ^° 



(A8) 



If the exact solution c(k, s) of Eq. I|A6(1 could be found, 
then the thermal-averaged velocity defined in Eq. © 



can be determined using standard limiting theorems for 
Laplace transforms as 



0c 



(v) = lim is — 
s^o dk 



(A9) 



k=0 



However, as Eq. (| A6|) is not in general exactly solv- 
able, we seek instead an approximate solution for c(k, s). 
A standard approach to approximating the solution of 
an integral equation is by iterating; i.e., first setting 
c(k, s) = P s (k) (neglecting the second term on the right 
hand of Eq. !A6Jl L then substituting this approximate 
solution into (|A6|) to obtain an updated solution, etc. 
After two iterations this procedure yields: 

c(k,s) = P s (k) - -^P s (k) J dpku(p)P(k p) 



(2tt)' 



P s (k) / dpku(p)P s (k p) 



dq(k-p) -u(q)P s (k-p-q) + 



(A10) 

with the dots signifying the further terms in the formal 
iteration series. Successive powers of the propagator ap- 
pear at each iteration, and therefore the approximation 
is assumed to be better for large F and large T. 

Consider now the velocity (v) in Eq. I|A9|I which results 
from using the approximation that terminates the series 
(|A10p by neglecting the dots. First, we differentiate c 
with respect to the component kj and evaluate at k = 0: 

dc 
dkj 



k=0 



P s (0) / dpu j (p)P s (-p) 



-(i P «(°) / dpuj(p)P s (-p) 

x J dq (-p) • u(q)P s (-p - q). (All) 

Following (|A9|) , we now multiply this equation by is 2 and 
examine the limit as s — > 0, and find that the first term 
of (|A11() yields Fj, i.e. the unmodified influence of the 
external force on the velocity. The second term of (|A11|) 
reduces in the limit to 

j^jdpujip) (lim S P s (-p)), 

which vanishes, since lim s _,o sP s (— p) = unless p = 
[by and "(°) = [ b y (EUl- To evaluate the limit 

of the third term of (|A11(I , we write it as 



lim , 

s^o { (2tt) 4 



s J dpUj(p)P s (-p) I dqp ■ u(q)P s (-p - q) 

(A12) 



and note that 

lim sP,.(-p)P s (-p - q) 

— lim 



o [s+p 2 T -ip ■ F] [s + (p + q) 2 T - i(p + q) • F] 
0, unless p = 0orp + q = 0. (A13) 
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Since Uj(p) is zero if p = [by Eq. l|A4(l ] . we conclude 
that the only non-zero contribution of (|A12|I occurs when 
q = p, giving the value 

(2^ / dp ^ (p) r^-i P -F p -" ( - p) - (A14) 

Together with the contribution from the first term of 
(|A11|) . we thus have the approximation 

<V> W F + (2^ / dP a(p) Tp»-q>-F P ' " ( " P) - 

(A15) 



Equation (JSJ) follows on replacing u from equation ilA3|l : 
the random potential case then requires a further ensem- 
ble average over the disorder to yield equation (fT2"|) . 
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